Correlations and phase structure of Ising models at complex temperature 
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We investigate the spin-spin correlation functions of Ising magnets at complex values of the 
temperature, T. For one-dimensional chain and ladder systems, we show the existence of a kind 
of helimagnetic order in the vicinity of contours where the leading two eigenvalues of the transfer 
matrix become equal in magnitude. We analyse the development of long-range order as the two- 
dimensional limit is approached, and find that there is rich structure in much of the complex-T plane. 
In particular, and contrary to the work of Fisher on this problem, the development of long-range 
order is actually associated with a proliferation of partition function zeros in a certain finite region 
of that plane containing the real-temperature magnetically ordered phase. The thermodynamic 
consequences of this are also discussed. 
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Introduction. There has recently been a resurgence of 
interest in the technique of using complex temperature 
as a probe of the physics of statistical mechanical models 
[IHl] , both in and out of equilibrium. The idea of extend- 
ing the equilibrium free energy to complex values of the 
temperature was first introduced by Fisher '^, following 
the related idea of using a complex magnetic field due 
to Lee and Yang [HITj. However, this has always been 
viewed as something of a technical trick, notwithstand- 
ing occasional statements in the literature such as "[e]ach 
region of holomorphy [of the free energy as a function of 
complex temperature] — if it comprises at least parts 
of the real axis — can be interpreted as a generalized 
'phase' of the system" [8]. 

Given that some of the abovementioned recent work 
[3] explicitly exploits the relationship between real time 
(relevant to dynamics) and imaginary temperature (as 
used by Fisher), it may be time to revisit the complex- 
temperature treatment of statistical mechanical models. 
In particular, a more solid interpretation of the physics 
at complex temperature should be provided. In what 
sense is each region of holomorphy a "generalized 'phase' 
of the system"? Is the qualification "if it comprises at 
least parts of the real axis" necessary? 

In this Letter, we address these questions by exam- 
ining the correlation functions at complex temperature, 
and the way in which long-range order in those corre- 
lation functions is related to the (more usually studied) 
distribution of zeros of the partition function. We do this 
in detail for a few particular cases: the one-dimensional 
Ising chain, quasi-one-dimensional Ising ladders, and the 
two-dimensional square-lattice Ising model. However, we 
believe that several of our conclusions are more generally 
applicable; we give arguments to this effect below. 

One- dimensional Ising chain. Let us begin by briefly 



analysing the one-dimensional Ising model at complex 
temperature. The Hamiltonian is 



H 



N 






(1) 



where the spin variable on each site aj = ±1, and we 
impose periodic boundary conditions so that cttv-i-i = ci- 
The partition function may be very easily obtained using 
a transfer-matrix method |9] : 



2^cosh^if + 2^sinh^if, 



(2) 



where K = J/ksT. The second term is often dropped 
because it is sub-dominant in the thermodynamic limit. 
However, it is cancellation between these two terms that 
provides the main structure of the model at complex tem- 
perature. In addition, as we shall see below, the sub- 
dominant transfer matrix eigenvalues determine the cor- 
relation length of the spin-spin correlator. 

The conventional approach is now to calculate the val- 
ues of K for which the partition function (pi) is zero; these 
correspond to singular points of the free energy, and thus 
are related to phase transitions. Clearly (pi) is zero when 



(tanhi^:)^ = -1 
or, taking the A^th root, 

tanhA: = e'^", 



(3) 



(4) 



where the angles On specify N equally-spaced points on 
the unit circle: 

^n = ^(" + ^), n = 0,l,2,...,A-l. (5) 



Hence the zeros of this partition function all lie on the 
unit circle in the complex tanhK plane. Since this vari- 
able is clearly the natural one when studying partition 
function zeros in Ising problems, we shall henceforth de- 
fine z = tanhK. We shall sometimes, for brevity, refer 
to this variable as the temperature, even though really it 
only encodes it. 

Non-analyticities in the free energy, i.e. phase transi- 
tions, can occur in general when we move from a region 
of the z-plane in which one transfer-matrix eigenvalue is 
dominant to a region in which a different one is. Hence, 
at a phase transition, we expect the two leading eigen- 
values of the transfer matrix to become degenerate. This 
criterion, however, is deeply related to the existence of ze- 
ros of the partition function. Indeed, there is a theorem 
due to Beraha, Kahane, and Weiss P^, which essentially 
states that (except for certain isolated points) functions 
of the form 



3 



(6) 



can have zeros only at points where |Ai| = IA2I, where 
the eigenvalues Xj are numbered in order of their magni- 
tudes with the largest first. Hence any phase transition 
must coincide with the crossing of a contour of partition 
function zeros. 

To characterise such transitions, we consider the latent 
heat. This can be evaluated as the difference in the in- 
ternal energy per particle between a point just 'inside' 
the unit circle and a point just 'outside' it. Continuing 
the usual thermodynamic formulas to complex values of 
the temperature, it is easy to show that 



AE 



Eout 



Ein 



= x: 



-1 oXin 



-x: 



-1 dXo 



dp 



(7) 



where Ain is the transfer matrix eigenvalue that is dom- 
inant inside the unit circle (2cosh_?C in this case), Aout 
is the one that is dominant outside (2sinhiir), and (3 = 
l/fcsT. Evaluating this latent heat at a temperature 



je 



on the unit circle, we obtain 



AE 



2iJ siiiO. 



(8) 



We see that, for real temperatures {6 = ot 6 = tt) 
there is no latent heat, and the transition is continuous; 
whereas for any other temperature the transition is first- 
order, with a purely imaginary latent heat. 

As well as determining the locations of the partition 
function zeros and of the phase transitions, eigenvalue de- 
generacies also determine the temperatures at which one 
sees long-range order in the correlation functions. This 
is because, in general, correlation functions are expressed 
in terms of transfer-matrix eigenvalues as 



(aoa,)=72(^) +73 (^) 



+ .. 



(9) 




FIG. 1. Contours in the complex tanh(J/fcsr) plane along 
which the two leading transfer-matrix eigenvalues are degen- 
erate for Ising ladders of various widths: (a) two legs; (b) 
three legs; (c) four legs. The dashed lines in panel (c) show, 
for comparison, the 'Fisher circles' [5] on which the partition 
function zeros of the isotropic two-dimensional Ising model 
lie. The solid circles on the real-temperature axis correspond 
to the ordering transitions of the isotropic two-dimensional 
Ising model: antiferromagnetic at 2 = 1 — y/2; ferromagnetic 
at z = ^2- 1. 



For this to be non-zero as j — ?► 00 , it must be the case that 
|'^2 1 = I Ai I , which is also the criterion for the possibility of 
a partition function zero (as per Beraha-Kahane- Weiss) 
and the existence of a non-analyticity in the free energy 
(see above). 

Notice, however, that (l9| is also sensitive to the rel- 
ative phase of the two largest eigenvalues. In the Ising 
chain, this gives rise to helimagnetic correlations when 
the temperature is complex. To see this, consider a tem- 
perature z = i?e'^ which lies near the unit circle. Explicit 
calculation gives 

C, = {cToa,) = 1 ^ JiN^^Ne • (10) 

For temperatures inside the unit circle, R is slightly less 
than one, so that for N ^ j we obtain Cj ~ We^^^ . 
Outside the unit circle, R is slightly greater than one, 
so that Cj ^ R^^e^'^^^. In both cases the correlations 
decay exponentially with distance, though the correlation 
length diverges as the temperature approaches the unit 
circle. But the correlations also have a spiral nature, 
and the handedness of this spiral switches as the circle is 
crossed. On the unit circle itself, the two spirals combine 



to form a density wave, 



C, ^ cos [Oj) . 



(11) 



Notice that this density wave evolves smoothly from uni- 
form long-range order at the ferromagnetic point {6 — 0) 
to Neel long-range order at the antiferromagnetic point 
(6* — tt) as the unit circle is traversed. Furthermore, the 
correlation functions are less sensitive to the boundary 
conditions of the problem than the contours of partition 
function zeros are, since the correlation functions sam- 
ple the sub-dominant eigenvalue of the transfer matrix 
even if that eigenvalue is accidentally absent from the 
partition function. 

One of the features of the treatment of Ising models 
at complex temperature is that the Z2 spins aj are ef- 
fectively promoted to C/(l) (i.e. XY-type) spins. This 
is not through a change in the nature of the partition 
sum, however, but through the fact that the statistical 
weights used in calculating the correlation functions are 
now complex. Hence it is strictly expectation values, not 
the spins themselves, that have their symmetry group 
enlarged. However, much of the intuition of spiral mag- 
netism of XY-spins at real temperature carries across to 
this complex-temperature Ising case. 

Ising ladders. It is natural to wonder which, if any, 
of these features survive into higher dimensions. To ad- 
dress this question, we now consider repeating the above 
analysis for the case of Ising ladders. Consider a ladder 
with Nl legs, and N :s> Nl sites on each leg. We take 
periodic boundary conditions in the 'long' direction, but 
open boundary conditions in the 'short' direction. For 
any number of legs greater than two, the problem already 
becomes analytically intractable, so most of our results 
in this section are numerical. We shall first present these, 
and then give an analytical derivation of some of the ob- 
served properties near the real-temperature line. 

Fig. [1] shows the contours along which |Ai| = IA2I for 
Ising ladders with two, three, and four legs. These are 
determined numerically, though of course for the Nl = 2 
case analytical calculation is also possible and agrees with 
the numerical results. In Fig. IlTc), as well as showing an 
expanded region of the complex-temperature plane, we 
have added the two circles along which Fisher |S] deter- 
mined the partition-function zeros of the isotropic two- 
dimensional Ising model to lie. 

Note the regular structure that appears around the 
points z = 1 and z = —1, which correspond respec- 
tively to the ferromagnetically and antiferromagnetically 
ordered T — states of the ladder. It seems from the 
figure that, for a ladder with Nl legs, precisely 2Nl de- 
generacy contours meet at each of these points, and that 
those contours are equally spaced in angle. This can 
in fact be proved as follows. Consider making a low- 
temperature expansion of the partition function around 
the point z = 1. The two dominant configurations of the 
system are the two ferromagnetic ground states, which 



are present at T = 0. As further terms in the expan- 
sion are included, these will dress the two ground states; 
but they will dress both of them in the same way, so 
the degeneracy between their contributions to the par- 
tition function will be maintained. The first excitation 
that breaks that degeneracy is a domain wall extending 
across the ladder in the short direction. Its effect may be 
described by the reduced transfer matrix 



w 



w 



Nl 
1 



(12) 



The entries on the diagonal are the contributions of the 
two ground states to the (renormalised) partition func- 
tion; each off-diagonal entry is the weight corresponding 
to a domain wall, where w = e"^''/'^^^ is the weight cor- 
responding to an unsatisfied bond. The eigenvalues of 
this matrix are 



A+ = 1 ± 



..Nr 



(13) 



and these are the dominant two eigenvalues of the full 
transfer matrix. They become degenerate when |A+| = 
|A_|; it is easily shown that this condition is satisfied 
when 



I — R exp 



ITT 



n=l,2. 



,2Nr 



(14) 

i.e. that 2Nl equally spaced degeneracy contours con- 
verge at the point z = 1. The analogous result for the 
z — —1 point may be established by very similar argu- 
ments. Note that we must have R -^ 1 for the low- 
temperature expansion to remain valid. 

We can also use this low- temperature expansion to ex- 
plore the latent heat associated with crossing each of 
these contours. Using (l7|, we obtain 



AE 

nn7, 



= ij(2-i?e*^)i?^^(-l) 



1 + i?2W^ ' 



(15) 



where 9 is one of the angles given in ( 14 ) . Notice that the 
latent heat vanishes as the d — 2 Ising limit {Nl — >■ 00) 
is taken. These observations continue to hold true even 
quite far from the real-temperature line: Fig. [2] shows 
the real and imaginary parts of the latent heat as the 
inner contour of Fig. fl|a) is crossed in an outward radial 
direction at angle 9 to the real line. 

We may also ask about the correlation functions in the 
ladder case. Fig. [3] shows the correlation function Cj for 
two spins on the same leg of the ladder with Nl = 2 and 
N = 100, separated in the leg direction by a distance j. 



We see that the density- wave behaviour of (11) persists. 



albeit with reduced amplitude, suggesting that our inter- 
pretation in terms of effective helimagnetism continues 
to apply in the ladder case. 

Approaching the two-dimensional Ising model. The 
above results give intriguing insights into the approach 



Re(AE) and Im(AE) 




\ 0,5 



FIG. 2. The real and imaginary parts of the latent heat per 
spin when the phase transition corresponding to the inner 
contour of Fig.fTfa) is crossed via a radial path in the complex 
temperature plane at an angle 9 to the real line. Notice that 
the latent heat is purely imaginary (as in the Ising chain case) 
until 9 reaches a critical angle corresponding to the kink in 
the contour in Fig. ITJa); after this a finite real part develops. 
We conjecture that this is a signature of a crossover from one- 
to two-dimensional behaviour. 
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FIG. 3. The spin-spin correlation function (uoaj) for two 
spins on the same leg of a two-leg ladder with A'^ = 100 spins 
on each leg. 



to the two-dimensional Ising model as N^ — > oo. Follow- 
ing the pattern of the preceding two sections, let us ask 
four questions: First, what happens to the eigenvalue- 
degeneracy contours in the Nl — t- oo limit? Second, 
what happens to the partition function zeros? Third, 
what happens to the latent heat? And fourth, what hap- 
pens to the correlation functions? 

We can confidently answer the first question near the 



zero-temperature point, since it follows from (14) that 



the region around z ~ 1 fills in entirely with such con- 
tours. Since the only significant lines in the complex 
temperature plane of the two-dimensional Ising model in 
the thermodynamic limit are the Fisher circles, it seems 
reasonable to conjecture that the entire region inside one 
Fisher circle but outside the other fills in likewise. This 



would imply that |Ai| — IA2I in this entire region in the 
Nl — >■ 00 limit. 

The answer to the second question involves a subtlety 
in the Beraha-Kahane- Weiss theorem: while degeneracy 
of eigenvalues is necessary for the existence of partition 
function zeros, it is not sufficient. In particular, if the 
condition |Ai| = jAa] is obeyed over a finite region in 
the z-plane, it does not lead to zeros of the partition 
function. But we do learn that the N^ —> 00 limit is 
somewhat singular in this regard. 

One might still ask, however, why Fisher's method [TT] 
based on Kasteleyn's solution of the dimer-covering prob- 
lem |12j does not see the extra partition function zeros 
in the case of a two-dimensional lattice of finite extent. 
There is nothing in the Fisher-Kasteleyn method that re- 
lies on the thermodynamic limit, so in principle all trans- 
fer matrix eigenvalues are there. We believe that the res- 
olution of this conundrum is related to the summation 
over different topological sectors that is required in the 
treatment of the two-dimensional Ising model. Fisher [S] 
in fact uses the partition function from only one of these 
sectors to determine the locations of the circles shown in 
Fig. [TJc). We conjecture that cancellations between the 
partition functions of the different sectors give the extra 
zeros that we find in the finite-size case. 



Regarding the third question, the formula ( 15 1 strongly 



suggests that, as the eigenvalue-degeneracy contours pro- 
liferate, the latent heat associated with each one is re- 
duced in magnitude (while retaining, in general, both 
real and imaginary parts). Indeed, we can see from (15) 
that even the sum of AE/NNl over all the transitions 
tends to zero as Nj^ — >■ 00, which implies that in the 
vicinity of the z = 1 point the thermodynamic behaviour 
becomes entirely regular. The only thermodynamic dis- 
continuity left is on the Fisher circles themselves, where 
the transition is presumably first-order. 

Finally, what of the correlation functions? There is 
nothing in the results obtained above to suggest that 
these revert to being purely ferromagnetic as the Nj^ -^ 
00 limit is taken; and on the other hand there is plenty of 
evidence [T3] to suggest that they retain a helimagnetic 
character. We conjecture that this is a general property 
of the correlation functions of Ising magnets at complex 
temperature: helimagnetism, with a pitch determined by 
the position of the temperature on the dense limit of the 
'co-ordinate grid' of contours seen in Fig. [T] 

This conjecture provides a natural framework in which 
to understand the instantaneous fragmentation of the re- 
gion into multiple separated phases as the lattice is made 
finite in one direction: it is a commensuration effect. Sim- 
ilarly, if one had a spiral magnet with a continuously 
evolving q-vector, finite-size effects would force the q- 
vector instead to make a sequence of transitions between 
the values permitted by the boundary conditions. 

Direct numerical or analytical evidence for these he- 
limagnetic correlations at complex temperature would 



be extremely useful. We also hope that these observa- 
tions may help in understanding the link between com- 
plex temperature and real-time dynamical problems. 

Conclusion. At the beginning of this Letter, we asked 
whether regions of holomorphy in the complex plane can 
properly be interpreted as generalized 'phases'. The 
answer for one- and quasi-one-dimensional systems is 
clearly no. Although these regions exist, long-range order 
is possible only on the degeneracy contours, which for un- 
frustrated systems are generally isolated lines. (For the 
case of frustrated models, see [13].) In two-dimensional 
models we have two different types of holomorphic re- 
gion in the complex-temperature plane: type I regions 
were holomorphic even in the quasi-one-dimensional ver- 
sion of the model, while type II regions became holo- 
morphic via the proliferation and coalescence of degen- 
eracy contours. In type II regions, we see helimagnetic 
order with a q-vector which evolves continuously as the 
region is traversed. The intersection of a type II region 
with the real axis corresponds to long-range order in the 
real-temperature model. Further work on this interesting 
classification of complex-temperature phases is desirable. 
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